UNIQUENESS AND NON-UNIQUENESS IN INVERSE RADIATIVE 

TRANSFER 



PLAMEN STEFANOV AND ALEXANDRU TAMASAN 

Abstract. We characterize the non-uniqueness in the inverse problem for the 
stationary transport model, in which the absorption a and the scattering coefficient 
k of the media are to be recovered from the albedo operator. We show that "gauge 
equivalent" pairs (a, k) yield the same albedo operator, and that we can recover 
uniquely the class of the gauge equivalent pairs. We apply this result to show 
unique determination of the media when the absorption a depends on the line of 
travel through each point while scattering k obeys a symmetry property. Previously 
known results concerned directional independent absorption a. 



1. Introduction 

This paper considers the problem of recovering the absorption and scattering prop- 
erties of a bounded, convex medium O C l n , n > 3 from the spatial-angular mea- 
surements of the density of particles at the boundary dfl. Provided that the particles 
interact with the medium but not with each other, the radiation transfer in the 
steady-state can be modeled by the transport equation 

(1) - 9 ■ Vu(x, 9) - a(x, 9)u(x, 9) + k(x, 9', 9)u(x, 9')d9' = 0, 

for x G Q and 9 G S 71 ' 1 ; see, e.g. 0, [22]. The function u(x, 9) represents the density 
of particles at x traveling in the direction 9, a(x, 9) is the absorption coefficient at x 
for particles traveling in the direction of 9, and k(x, 9', 9) is the scattering coefficient 
(or the collision kernel) which accounts for particles from an arbitrary direction 9' 
which scatter in the direction of travel 9. 

The medium is probed with the given radiation 

(2) u\ r _ = /_ 
and the exiting radiation 

(3) u\ r+ =: A[fJ\ 
is detected; where 

(4) r± = {{x, 9) G dQ x S^ 1 : ±6 ■ n(x) > 0} 

with n(x) denoting the outer unit normal at a boundary point x. The equation ([3]) 
defines the albedo operator A which takes the incoming flux /_ to the outgoing flux 
u\r + at the boundary. 
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In general, the boundary value problem ([T|) and (j2J) may not be uniquely solvable 
but it has a unique solution for generic (a,k); see [2E] and Proposition 13.11 below. 
Unique solvability can also be obtained under some physically relevant subcritical 
conditions like (jEJ) or ([7J below. We assume that (a, k) is such that the direct problem 
(PQ) and is well posed. 

The inverse boundary value problem of radiative transfer is to recover the absorp- 
tion a(x, 9) and the scattering kernel k(x, 9', 9) from knowledge of the albedo operator 
A. One could also study the time-dependent version of (TO , and then the kernel of A 
contains one more variable that gives us more information. This problem has been 
solved under some restrictive assumptions (e.g. A; of a special type or independent of 
a variable) in [U [2J [31 [SJ [TH [13 Q21 [271 EH] • In three or higher dimensions, uniqueness 
and reconstruction results for general k and a = a(x) were established in [TT]. The 
approach there is based on the study of the singularities of the fundamental solution 
of ([1]) (see also [8]), and the singularities of the Schwartz kernel of A. Stability es- 
timates for k of special type were established in [2~T1 [29] ; and recently, for general 
k, in j6]. Uniqueness and reconstruction results in a Riemannian geometry setting, 
including recovery of a simple metric, were established in [19]. Similar results for 
the time-dependent model were established in [10], and in [13] for the Riemannian 
case. In planar domains the work in [25] shows stable determination of the isotropic 
absorption and small scattering, and extensions to simple Riemannian geometry is 
given in [2tJ]. Also in two dimensional domains we point out that the recovery of 
k is only known under smallness conditions which are more restrictive than what 
is needed to solve the direct problem; e.g. more restrictive than ([51) or (17)) below. 
On the other hand, in the time-dependent case, the extra variable allows us to treat 
the planar case without such restrictions, see [10]. We also mention here the recent 
works [7J [HI [15], in which the coefficients are recovered from angularly averaged 
measurements rather than from knowledge of the whole the albedo operator A. 

The above mentioned results concern media with directionally independent absorp- 
tion property; for transport with variable velocity v, which now belongs to an open 
subset of M. n , the absorption may depend on the speed a = a(x, \v\) . 

In general, in media with an anisotropic absorption, the albedo operator does 
not determine the coefficients uniquely. For example, if k = the obstruction to 
unique determination can be readily seen. The most one can recover from the albedo 
operator A are the integrals La(x + t9,9)dt, see [TT]. In [6], they are shown to be 
stably recovered independently of k. In other words, for each fixed direction 9 G S n ~ 1 , 
we know the integral of the map a(-,9) over the parallel lines in the direction of 9. 
This is insufficient data, since one can smoothly change the x-variable in the direction 
of 9 while preserving the integral; see also the theorem below. 

The non-uniqueness described above and the need to assume that a is isotropic 
left the theory of the inverse radiative problem in a somewhat unsatisfactory state. 
It was not clear whether the uniqueness failed if k ^ 0, neither it was known what 
information about (a, k) can still be recovered. The purpose of this work is to fill 
this gap. We show that a certain family of "gauge transformations" of (a, k) does 
not change the albedo operator A; and that, given A, one can recover uniquely the 
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class of gauge equivalent pairs (a,k). The recovery of the gauge equivalent class is 
explicit, as the recovery of a = a(x) and k in [TTJ is explicit. 

2. Main Results 
The pair (a, k) is assumed to satisfy the admissibility condition 

(5) sup (\a{x,6)\+ [ \k{x,9,9')\d9') < oo. 

Let T be the operator defined by the l.h.s. of ([TJ). For (a, k) admissible, the second 
and the third terms of T are bounded operators on L x (f2 x >S n_1 ), while the first 
term is unbounded. We view T as an unbounded operator on L x (f2 x S"™ -1 ) with the 
domain 

D(T) = {f e L\n x S n ~ l ); 6-VfE L\n x S^ 1 ), f\ r _ = 0}, 

see also [TTJ. Since the direct problem (pQ) and (j2J) can always be reduced to a non- 
homogeneous problem with a homogeneous boundary condition, invertibility of T 
implies well-posedness. We say that the direct problem is well posed, if T~ l exists as 
a bounded operator. 

As an example, we have the following two subcritical conditions that yield well- 
posedness; see, e.g., [HI HJJ CE21 ED [22] and Proposition 13.11 below. Either 



(6) sup 



t(x,6) / k(x,6,6')d6 f 



< 1, 



where t(x,9) is the total free path of a particle at (x,9), see the beginning of the 
next section, or 



(7) a(x, 



[ k(x, 9, 9')d9' > 0, a.e. Q x S" -1 . 



both a and A; outside Q to 

We start with a simple observation, that seems to be new. We will show that 
there is non-uniqueness even if k ^ 0, too. Let <p(x, 9) > be such that = 1 on 
T := dQ x S' n ~ 1 . Set a = a — 9 ■ V x log 0. Then we can rewrite ([1]) as 

-9 -X7 x u- <p~\9 -V x <p)u-au+ [ k(-, 9', 9)u(-, 9')d9' = 0. 

Multiply by to get 

-9 ■ V x ((pu) -acfm + cj) [ k{-, 9', >(-, = 0. 

The function u — cjm thus solves (pQ) with (a, k) replaced by (a, fc), where those two 
pairs are related by the "gauge transformation" 

(8) a = a-9- V x .log0, k(x,9',9) = ^^j-k(x,9' ,9). 

Since ^ = 1 on T, the boundary data do not change. Therefore, if A is the albedo 
operator corresponding to the pair (a, k), then A = A. 
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Our main result is that this is the only obstruction to non-uniqueness. 

Theorem 2.1. Let (a,k) and (a,k) be two admissible pairs for which the direct 
problem is well-posed, and let A and A be the corresponding albedo operators. Then 
A = A if and only if there exists a positive <fi 6 L°°(Q x S™' 1 ) with 9 ■ V x (p(x, 9) € 
L°°(tt x S' 11 - 1 ) and$=\ for x E dtt, so that © hold. 

The proof of the theorem is based on the analysis of the singularities of the Schwartz 
kernel of A, as in [TT] . 

Theorem 12.11 allows us to obtain a few new uniqueness results under additional 
conditions. One of them concerns the case where we have anisotropic media with 
absorption a(x, 9) that depends on the line of travel through each point (but not on 
the direction): 

(9) a(x, 9) = a(x, -9), xgO, 9 e S n -\ 

and a scattering coefficient k > satisfying the following symmetry condition 

(10) k(x, 9, 9') = k(x, 9', 9), xett, 0,9' € S"' 1 . 

Corollary 2.2. Let (a,k), (a,k) be two admissible and subcritical pairs which yield 
the same albedo operator. Assume that k and k > satisfy f fTOj) . 

(i) Then k = k and a = a + 9 ■ Vv(x) for some function v(x) vanishing on 
dVt. In particular, one can recover the total absorption at a.e. x, i.e., J a(x,9)d9 = 
Ja(x,9)d9. 

(ii) If, in addition, a and a satisfy then a = a. 

Note that any two pairs as in (i) yield the same albedo operator, so this is the most 
we can say in this case. The symmetry assumption (TTDT) occurs naturally in some 
models of Optical Tomography, where the scattering of light in a tissue depends on 
the angle between the two directions: k(x, 9, 9') = k(x, 9-9'). 

One can formulate and prove similar results in the case where the velocity belongs 
to an open subspace of R n , i.e., the speed can change, as in [TT]. We restrict ourselves 
to the fixed speed case (\9\ = 1) for the sake of simplicity of the exposition. Also, the 
fixed speed model is the one that is most often discussed in the literature. 

The paper is organized as follows. Section [3] recalls some results from [TT] that 
we use later. In Section H] we prove Theorem 12.11 and its Corollary 12.21 In section 
[5] we give the reconstruction formulae for continuous a and k in the symmetric case 
covered by Corollary 12.21 Section [Tj] contains concluding remarks. 

3. Preliminaries 

In this section we recall some results from [TT] recast to the one-speed velocities 
and introduce notations. 

Let t±(x,9) be the travel time it takes a particle at x to reach the boundary 
while moving in the direction of ±9 and define t(x,9) = r_(x,6*) + t + (x,9). Since 
we work with unit-speed velocities, note that t(x,9) < diam(f2). Let d^{x,9) = 
\n(x) ■ 9\dfi(x)d9 where d/i(x) is the induced Lebesgue measure on the boundary and 
d9 is the normalized measure on the sphere. Also let 5{ x }(x') represent the delta 
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distribution with respect to the boundary measure dfj,(x') supported at x G dQ and 
5{g}(9') represents the delta distribution with respect to d6 centered at 9 G >S n_1 . 

Also, for (x,9) G Q x S n ~ 1 , let x\ denote the exiting point if traveling from x in 
the ^-direction and Xg, be the entrance point at the boundary to reach the inside 
point x by traveling in the ^'-direction, i.e. 

(11) x^ := x + t + (x, 9) and Xq, := x — r_(x, 6')6' . 

For the proposition below, we introduce the class of regular scattering kernels 
k G C(Q, L^iS^^L^S 71 - 1 ))). Then the map / k(x,9' ,9)4>{9')d9' is bounded on 
L 1 (S' n ~ 1 ) continuously depending on x. Our notion of regular k is stronger than that 
in [21], and, in particular, it allows us to use the results in there. 

Proposition 3.1. The direct problem is well-posed, i.e., T has a bounded inverse on 
L l (Q x S 17 " 1 ), if either or ([7]) holds. Moreover, the direct problem is well-posed 
for an open dense set of(a,k) G L°°(Q x S^ 1 ) x C(fi, L°° (S^ 1 , L 1 (S^ 1 ))) . 

Proof. We will first discuss the well-posedness under the conditions (EJ) or ([7]) . Assume 
§6§ first. By [TJ1 Proposition 2.3], and since r is bounded, we get that T~ l is bounded. 
The subcritical case (J7J) is covered in P, Section 2]. 

The generic statement is proven in L 2 (Q x S'™ -1 ) for C 2 coefficients in [26]. In the 
L 1 spaces under consideration, we proceed in a similar way. Let K be the integral 
operator in ([1]) and T\ = T — K. Then KT^ l K is weakly compact in L l {Vt x S™" 1 ), 
see [2T]. Therefore, (KT^ 1 ) 2 is weakly compact, and its square is compact. For 
a fixed a, consider the family AA;, where A is a real parameter. By the analytic 
Fredholm alternative in Banach spaces [23], A i— ► (/— (AA'Tf 1 ) 4 ) -1 is a meromorphic 
family. This implies that A t— > (I — XKT^ 1 )" 1 is also meromorphic, and thus T^ 1 = 
(/ — AA'T]" 1 ) -1 ^" 1 exists for all but a discrete set of A's. This shows that there is a 
dense set of pairs yielding a well-posed problem. The fact that this set is also open 
follows from a perturbation argument around each (a, k), for which T _1 is bounded; 
thus (/ — AATf 1 ) -1 corresponding to nearby pairs exists. □ 

Note that one can set a(x,9) = ao(x,9) + J k(x,9,9')d9', where the integral rep- 
resents the attenuation due to the change of direction, while is the absorption. 
Then one can prove in the same way that the direct problem is well posed for generic 
(ao, k), and moreover, for any a , this is true for generic &'s. 

Proposition 3.2. Assume that the direct problem is well-posed. Then the albedo 
operator A : L 1 (r„,(i^) — > L 1 (r + ,(i^) is bounded and its Schwartz kernel a is given 
by a = a± + a2 + «3, where 



ai(x,9,x',9') =ef° 
ct2(x, 9, x , 9') = I 



T_(x,e) 



x - te)dt 5 {x - } {x')5 {e} {9') 



(12) 



T-(X,6) 



e - J£a(x-tO,0)dt e - f t 



T_(x~se,e') 



(13) 



a{x-s0-t9' ,6')dt 



Jo 



(14) 



x k(x - s9, 9', 9)5 {x _ s e_ T _ {x _ s e t ei)o>}{x')ds 
n{x') ■ ff^ats GL°°(r_; L 1 (T + , d£)). 
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This proposition is formulated in [TT] under the assumption that the system is 
subcritical, i.e., either or (j7|) holds. We remark that those conditions are only 
used in the proof in the analysis of 0:3, to guarantee that T _1 exists in L ; something 
that we assume here. 

Let <j) G Cq°(B(0; 1)) with < < 1 and 0=1 near the origin be a cut-off 
function. Given e > we define for x, x' G W 1 and 9, 9' G S n ~ l 

(15) <j> e (x,e,x',e') = <t> f 

Proposition 3.3. Assume that the direct problem is well posed. Then the limit below 
holds in L 1 (T + , d£) 

(16) lim / a(x,e,x',e f )(P e (x,9,x',e')dfi(x')de = e-^ a{x+t8 ' e)dt > a.e. T+; 

For linearly independent 6, 9' G S 71-1 , let ng^^x) denote the projection of x onto the 
plane through the origin spanned by 9 and 9'. Let 9' ± be a unit vector in span{#, 9'} 
orthogonal to 9': 9' ± ■ 9' = 0. Let cp G Cg°(-1, 1) with < p < 1 and / R ip{t)dt = 1 be 
a cut-off function in R and be the cut-off function in IR n introduced earlier. Define 
the test function 

(17) <p ej5 (x;9,9) = -Lp 



With xj" and given by ( TTTT) and (j) e> g given by ( 1TT1) . we define in fl x S 11 ^ 1 x S 1 ™" 1 



(18) I e> s(x,9' ,9) :— / a(x^,9,x',9')(j) ei s(x' — x e ,,9,9')dfji(x'). 

JdQ, 

Proposition 3.4. Assume that the direct problem is well posed. Then the following 
limit holds in L x (fi; L\ oc {F)) 

(19) IUx,ff,0) := \im\iml e5 (x,9',9) = e - /"„ «(*+*W-J " «(*+*^ 

where P := {(9,9') G S^ 1 x S™" 1 : 0' ^ ±9}. 

4. Proof of Theorem 12.11 and Corollary 12.21 
Proof of Theorem \2.1[ From Proposition 13. 3^ we can recover 

/oo 
a(x + s9,9)ds, (x,9) G tt x S n -\ 
-oo 

In particular, the integrals above for / = a — a vanish: 

/oo 
f(x + s9,9)ds = 0, V(z, 9) G x S n - X . 
■00 

The kernel of the linear transform (1201) is easy to describe. Set 



(22) v(x,6) 



/ f(x + s9,9)ds, (x, 9) E £1 x S" 2-1 . 
J —00 
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Then v E L°°(fi x S™" 1 ) with 9 ■ V x v = f E L°°(tt x S"- 1 ) and v = on T_. From 
(l2~Tj) we get also that i> = on T + , therefore v (x, 9) = for x E dfl and a.e. 9 E S n ~ 1 . 
Set 

(23) (j){x, 9) = e v{x ' 9) 

to get a — a = 9 ■ V x log <fi as claimed. This shows the first part of ([H]). 
Now, from Proposition 13.41 we get 

Using the first equality in (jHJ) that we already proved, after a simple calculation, we 
derive 

(24) k(x,9',9) = ^j- ) k(x,9' 1 9). 

□ 

Proof of Corollary ' \2.2l Next we use the characterization above to prove uniqueness 
in the symmetric case. By swapping 9 and 9' in (l24"j) we get 

For the second and fourth equality we used ([TO]) , while for the third one we used 
(l24"j) . Since /c does not vanish, we conclude that <j>{x, 9') = <fr(x, 9) for all 9, 9' E iS™ -1 . 
Therefore <j) = 4>(x) is independent of 9 and applying (l24"j) again we conclude that 
fc = A. 

So far we showed that a(x, 9) — a(x, 9) = 9 ■ V x (p(x). If the symmetry relation (Q 
holds then a(x, 9) — a(x, 9) = —9 ■ V x (p(x). Therefore a — a = 0. □ 

5. Reconstruction formulas in the symmetric case 

We showed in the previous section that under the symmetry hypotheses f fTUj) and 
Qj there is uniqueness for a, k. The proof is constructive but it still leads to the 
problem of recovering a(x,9) first, up to 9 ■ Vlog</>, from its under-determined X- 
ray transform, see the l.h.s. of (jHJ). In this section, we show that under the same 
assumptions, there is an explicit reconstruction of a different, simpler type. Namely, 
we recover the integrals in (1211) first, that measure the total absorption along straight 
lines, but we are not trying to determine a from them. Then by ([TBI) we recover 

see ( 1281) below. We do not know the attenuation term above because we know the 
attenuation along straight but not broken lines. We can however swap 9 and 9' and 
multiply the results. Then we get k 2 multiplied by the attenuation along the two 
lines through x parallel to 9, and 9', respectively, and we know that attenuation. 
This recovers k, and then we recover a. 
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We need to strengthen the regularity assumptions on the coefficients to 

(25) aeC(Qx S^ 1 ), k G C(Q x S^ 1 x S^ 1 ). 

Firstly, we extend the limit in Proposition 13.31 valid for maps in r + to maps in 
Q x S"' 1 . For (x, 9) G Q x S" 1 " 1 and e > 0, let us denote 



(26) 



J e (x,9) : = I a(xg ,9,x' ,8')<f) e (x,9,x' ,8')dfi(x')d8' 



where x^ is given in (fTTl) . We remark here that J e is constant with x varying in the 
direction of 9, i.e. J e (x + t9, 9) = J e (x, 9) for all t. 

Lemma 5.1. The limit below holds in L 1 ^ x 5" 1-1 ; t(x, 9)~ 1 dxd9) 
(27) J (x, 9) := lim J e (x, 9) = exp (- J a(x + t9, 9)dtj > a.e. 0x5' 



n-l 



Proof. We use the change of variable formula 



f(x, 9)dxd9 



nxs™ 



r+ Jo 



f(x+ - s9)dsd£(x + ,9) 



for f(x, 9) = J t (x, 9) — exp ^ j 00 ^ a(x + t9, 9)dt s j r(x, 9) 1 to estimate 



/ 



r+ Jo 



J e (x,9) - exp / a(x + t9,9)dt 



-(*+,*) 



T(x,9)- 1 dxd9 



J e (x+ - s9,9) - e {l-^{* + sO + te m ) 



J t (x + , 9) — exp 



r(x+ - s6>,6<) 
a(x + + t9,9)dt 



dsd£(x + ,9) 



The last identity holds due to the fact that t(x + — s9, 9) = t(x + , #) = r_(x + , 6*) and 
that the integrand is constant in s. An application of Proposition 13.31 finishes the 
proof. □ 

Following from Proposition 13.41 for 9, 9') G f2 x I" we get that 

(28) e -;° 00 «(*+^,«')*-; 00 «(*+^)*jfe( a;j ^^) = I 0fi ( x ,9',9), 

(29) e -/° oo a( a ,+t9,fl)*-/ oo a(x+W^)* A .( a . j ^ 0/) = j 0j0 ( a . j ^ ^ 

Using the symmetry assumption ffTUl) on k, by multiplication of (128]) and (1291) we 
obtain 

(30) e -n o «(*+« , .» , )*-r o o^+^)*fc2(x, 0', 0) = I ,o(x, 9', e)I 0lQ (x, 9, 9'). 
Now using Lemma [5.11 we recover 

I 0>0 (x,6',8)I 0>0 (x,9,6') 



(31) 



k(x,e',9) 



J {x,8)J {x,8>) 
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Using the formula (J3TJ) in ([29]) we get for a.e. (x, 9, 9') G Q x I" that 
(32) 

f a(x + t9,9)dt + r a (x + te',6')dt= Ibgf /;;f?f'f f/ .A 

7-oo Jo 2 \I 0fi (x,9,9')J {x,9)J {x,9>) J 

The continuity assumption (1251) imply that the identity above extends pointwise in 
f2 x S" 1 " 1 x S" 1 " 1 . By looking at backscattering 9' = —9 and making one change of 
variables t «-> — £ we obtain 

(33) 

£ a(I + *, ^< + jl a(x + *, = i log ( /M(Zifli ^ ( ;^,- fl ) ) ■ 

The left hand side is differentiable in x in the direction of 9, hence so is the right 
hand side. By differentiating ( 1331) in the direction of 9 we obtain 

(34) a (x ,e ) + a (x , -0) = \ 9 . V, log ( w , _,, ) ■ 

Since the absorption depends on the line through x, we have a(x, 9) = a(x, —9) and 
thus it can be recovered from the formula: 

(n\ ( a\ l fi T? i f Io,o(x, -9,9) \ 

35 a ( x > g ) = T^- V^log — ^ . 

4 \Iofi{x,9,-9)Jo{x,9)J {x,-9)J 

6. Remarks 

6.1. The isotropic absorption case. The previously known uniqueness result [11] 
for isotropic absorption follow from Theorem 12. II combined with the injectivity of the 
X-ray transform. If / = a — a depends on the position only, then (121!) implies / = 0. 
From the definition (T23|) we get = 1, which by (JHJ) yields k = k. 

6.2. Other conditions that imply uniqueness. Assuming ffTUl) . we can only re- 
cover a up to 9 ■ Vi>(x) as above. For unique recovery, the condition d^J) suffices 
but it can be replaced by something weaker. For example, we may require that a is 
orthogonal to all such functions w.r.t. some measure djj,{9). This is equivalent to 

/ 9 ■ V x a{x, e)df Jl {6) = 0. 

JS"- 1 

Note that the symmetry assumption (Q implies such a condition if dfi(9) is even. 
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